Role of conformational entropy in force-induced bio-polymer unfolding 
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A statistical mechanical description of flexible and semi-flexible polymer chains in a poor solvent is 
developed in the constant force and constant distance ensembles. We predict the existence of many 
intermediate states at low temperatures stabilized by the force. A unified response to pulling and 
compressing forces has been obtained in the constant distance ensemble. We show the signature of 
a cross-over length which increases linearly with the chain length. Below this cross-over length, the 
critical force of unfolding decreases with temperature, while above, it increases with temperature. 
For stiff chains, we report for the first time "saw-tooth" like behavior in the force-extension curves 
which has been seen earlier in the case of protein unfolding. 
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During the past decade force has been used as a ther- 
modynamic variable to understand molecular interac- 
tions and their role in the structure of bio-molecules 
[E 0j- By exerting a force in the pN range one can 
experimentally study the elastic, mechanical, structural 
and functional properties of bio-molecules [3[. The de- 
pendence of force on concentration, pH of the solvent, 
loading rate and temperature provides basic understand- 
ing of the interactions 0, H, @]. Many biological reac- 
tions involve large conformational changes which provide 
well defined mechanical reaction co-ordinates, e.g. the 
end-to-end distance of a polymer, that can be used to 
follow the progress of the reaction Q. Such processes 
have been modeled by a simple two state model Q . The 
applied force "tilts" the free energy surface along the re- 
action co-ordinate by an amount linearly dependent on 
the end-to-end distance. The kind of transitions induced 
by the applied force are the folding-unfolding transition 
of proteins 111 , the stretching and unzipping transition of 
dsDNA 0, [1] or the ball-string transition of a polymer 
[Tfj| . From polymer theory we know that a polymer chain 
will, depending on temperature, be in either a collapsed 
state or a swollen state [ll|. The end-to-end distance 
(R) scales as N v , where N is the chain length and v is 
the end-to-end distance exponent. In the collapsed state 
(low temperatures) v — 1/d, while at high temperatures 
v is given quite accurately by the Flory approximation 
v = 3/(d + 2) [11] (actually this formula is exact for 
d = 1, 2 and 4). ft should be noted that by varying 
temperature alone, a polymer can not acquire the con- 
formation of a stretched state where v = 1. Hence force 
not only "tilts" the free energy surface but also induces 
a new "stretched state" which is otherwise not accessi- 
ble. Moreover, recent experiments suggest that there are 
many intermediate states involved which are crucial to 
the understanding of unfolding experiments and lie be- 
yond the scope of two state models [l(| EH [IH . 

The non-equilibrium thermodynamics of small systems 
has mostly been studied in the "constant force ensemble" 



(CFE) where the control parameter is the average exten- 
sion. Most applications of atomic force microscopy apply 
the force using a linear ramp protocol with a very small 
velocity. Such systems may be considered to be in quasi- 
static equilibrium and the appropriate ensemble is the 
"constant distance ensemble" (CDE). In the thermody- 
namic limit both ensembles are expected to give similar 
results [14]. However, single molecule experiments study 
systems of finite size and the results may depend on the 
ensemble 15j. Apart from this, the physical constraints 
imposed by experimental setups have not been fully taken 
into account. For example, in atomic force microscopes, 
receptor and ligand molecules are attached to a substrate 
and a transducer, respectively. The loss of entropy due 
to the confinement has been ignored in most models. 

Theoretically these transitions have been studied us- 
ing simple models such as freely jointed chain (FJC) 
or worm like chain (WLC) models [H B3|. The WLC 
model has been used to study the force-extension curves 
of bio-molecules. However, this model ignores excluded 
volume effects and attractive interactions between chain 
segments and is thus only well suited for modeling the 
stretching of polymers in a good solvent [ItJ • For a poly- 
mer in a poor solvent, the force-extension curve shows 
the existence of a plateau region at a well defined force 
floL [l2| . Experimental observations of polymers in poor 
solvents and an improved theoretical understanding of 
semi-flexible polymers have given rise to theories of glob- 
ules with well defined internal structure [H, [l9| . These 
theories have potential applications in the study and un- 
derstanding of the basic mechanisms of protein folding. 

The purpose of this communication is firstly to pro- 
vide a complete phase diagram using exact results of fi- 
nite chains. However, we focus our studies on the be- 
havior of the force-extension curves at low temperatures 
where the usual Monte Carlo approach and the theoret- 
ical models discussed above fail. These studies have di- 
rect relevance to phenomena seen in experiments with 
single biomolecules. We report the existence of interme- 
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FIG. 1: An ISAW on the square lattice with one end attached 
to a surface and subject to a pulling force on the other end. 

diate states stabilized by the force at low temperatures. 
We also show the signature of a crossover length, below 
which the transition temperature decreases with force, 
while above the crossover length the transition tempera- 
ture increases with force. When bending rigidity is taken 
into account, we predict "saw-tooth" like oscillations in 
the regime of large polymer stiffness, as observed in pro- 
tein unfolding experiments. 

To obtain exact results for small chains we model them 
as interacting self-avoiding walks (ISAWs) on the square 
lattice [2(| as shown in Fig.Q] Interactions are introduced 
between non-bonded nearest neighbor monomers. In our 
model one end of the polymer is attached to an impene- 
trable neutral surface (there are no interactions with this 
surface) while the polymer is being pulled from the other 
end with a force acting along the x-axis to model the 
experimental setup. We introduce Boltzmann weights 
u> = exp(— e/k^T) and u = exp(— F/k^T) conjugate to 
the nearest neighbor interactions and force, respectively, 
where e is the interaction energy, fee is Boltzmann's con- 
stant, T the temperature and F the applied force. In the 
rest of this study we set e — — 1 and fee = 1- We study 
the finite-length partition functions 

Z N (F,T) = Yy m u x = ^2C(N,m,x)uj m u x , (1) 

all walks m.x 

where C(N,m,x) is the number of ISAWs of length 
N having m nearest neighbor contacts and whose 
end-points are a distance x — xn — xq apart. 
The partition functions of the CFE, Zn(F,T), and 
CDE, Z N {x,T) = £ m C(iV,m,2:)w m , are related by 
Zn(F,T) = J2 x Zn(x,T)u x . The free energies are eval- 
uated from the partition functions 

GO) = -T\ogZ N {x) and G(F) = -T log Z N (F). (2) 

Here (x) = dG Qp ^ and (F) = dG d ^ are the control pa- 
rameters of the CFE and CDE, respectively. 

We enumerate all possible conformations of the ISAW 
by exact enumeration techniques. The major advantage 
of this approach is that the complete finite-length par- 
tition functions can be analyzed exactly. Furthermore 
scaling corrections can be taken into account by suitable 
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FIG. 2: The phase diagram for flexible chains. 



extrapolation schemes enabling us to obtain accurate es- 
timates in the thermodynamic (infinite length) limit [2l| . 
To achieve a similar degree of accuracy using Monte Carlo 
simulations one typically has to use chains at least two 
orders of magnitude longer than in the exact enumera- 
tions |22j]. The greatest challenge facing exact enumer- 
ations is to increase the chain length. Until now most 
exact results for models of small proteins were confined 
to chain lengths of 30 or so [23|, Hj] • Here the number of 
ISAWs was calculated using transfer matrix techniques 
[25l | . Combined with parallel processing, we were able to 
almost double the chain length. To be precise we calcu- 
late the partition functions up to chain length 55. 

In Fig. [2j we show the force-temperature phase di- 
agram for flexible chains. The qualitative phase dia- 
gram remains largely the same [241 ] . The observed re- 
entrant behaviour can be explained by a ground state 
with non-zero entropy [241. l26j|. However, the extension 
up to length 55, allowed us to obtain a new transition 
line from the extended state to the stretched state which 
is solely induced by the applied force. In contrast to 
the lower phase boundary (collapse transition), where 
the force decreases with temperature, the upper phase 
boundary (stretching transition) shows that the force in- 
creases with the temperature. 

First we study the model in the CFE. In Fig. [3^, we 
plot the average scaled elongation with force for different 
chain lengths at low temperature. Experimentally sev- 
eral transitions were found in the force-extension curves 
corresponding to many intermediate states [9I ITol. [iH IT3I] . 
This phenomenon is clearly confirmed by our study. It 
has been argued that in the limit of infinite chain length, 
the intermediate states should vanish and there will be 
an abrupt transition between a folded and fully extended 
state [231 ]. Evidence to this effect can also be seen in 
Fig. [3ji where we note that the plateaus at an exten- 
sion around 0.2 tends to increase with N while the other 
plateaus tend to shrink with N (this is particularly so for 
the plateau around 0.5 corresponding to a simple zig-zag 
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FIG. 3: The average scaled elongation {x)/N vs F at T 
0.01 for various lengths (a) and temperatures at length N 
55 (b). 



FIG. 4: Plot of the average force (F) vs the elongation x 
at various temperatures T for N — 55 (a) and the cross-over 
length vs N (b). 



pattern of the chain). Note however, that as we change 
the chain length from 25 to 55, we find more and more of 
these intermediate states. This has also been observed in 
recent experiments [H, [1(3] where the globule deforms into 
an ellipse and then into a cylinder. At a critical exten- 
sion the polymer undergoes a sharp first order transition 
into a "ball string" conformation [{j, [l(| ■ This shows that 
finite size effects are crucial in all the single molecule ex- 
periments and can be seen even for long chains [1 31 ] . A 
simple theoretical argument for the observed behavior 
is that at low temperature, where the entropy S (per 
monomer) of the chain is quite low, the dominant contri- 
bution to the free energy 



G = Ne — a(N, F)e — NTS 



(3) 



is the non-bonded nearest-neighbor interaction Ne. The 
second term is a surface correction and it vanishes in the 
thermodynamic limit. However for finite N, the system 
has many degenerate states depending upon the shape 
of the globule. This leads to a surface correction term 
a(N,F) which is a function of N and F. If F = 
the shape of the globule is like a square and the sur- 
face correction term a(N, 0) will be minimized and equal 
to 2>/N. In the CFE, there is a force induced additional 
contribution proportional to the extension of the glob- 
ule which along with a(N, F) stabilizes the intermediate 
states. When the temperature increases the multi-step 
character of the force-extension curve is washed out due 
to increased contributions from entropy. This effect can 
be clearly seen in Fig. [5Jd where we have plotted force- 
extension curves at different T. 

Next we study the model in the CDE. The force- 
extension curve shown in the insert of Fig. \Ql has inter- 
esting features. It shows that when the distance between 
the first and the last monomer (where force is applied) 
is less than the average size of the coil (without force), 
one needs a compressing force instead of a pulling force. 
The qualitative behavior is similar to one observed in 
experiments [13] and computer simulations [28|. Since 
most models do not include confinement in their descrip- 
tion, such behavior could not be predicted. In Fig. 3^,, 
we show the response of the force when the elongation 



exceeds the average size of the polymer. The flat por- 
tion of the curve gives the average force needed to unfold 
the chain. Such plateaus have been seen in experiments 
Ed, mill. From Fig. Hi one can also see that the force 
required to obtain a given extension initially decreases 
with temperature. But beyond a certain extension (close 
to 30 in this case) the required force increases with tem- 
perature. We note that the curves cross each other at 
a 'critical' extension for any temperature (below the 9- 
point). We identify this as a cross-over point. In Fig. [Dd 
we plot the position of the cross-over point L cross as a 
function of the length N of the polymer chain. We see 
that the crossover extension increases linearly with the 
chain length. This shows that above this point the chain 
acquires the conformation of the stretched state. The in- 
crease in force with temperature generates a tension in 
the chain sufficient to overcome the entropic effect. Since 
the contribution to the free energy from this term is TS 
(S being the entropy), more force is needed at higher T as 
seen in experiments. Our exact analysis for finite chain 
length shows that applying a force at first favors taking 
the polymer from the folded state to the unfolded state. 
However, rupture or second unfolding occurs when the 
tethered or unfolded chain attains the stretched state, 
and one requires more force at higher temperature. 

We model semi-flexible polymers by associating a pos- 
itive energy ej, with each turn or bend of the walk 24j • 
The corresponding Boltzmann weight is loi, — exp(— bet), 
where b is the number of bends in the ISAW. We again 
enumerate all walks, but because of the additional pa- 
rameter uib, we were restricted to 45 steps. For a semi- 
flexible polymer chain, a stretched state may be favored 
by increasing the stiffness. The phase diagram for semi- 
flexible chains is now well established. It has three states 
namely (i) an open coil state at high temperature, (ii) a 
molten globule at low temperature and low stiffness and 
(iii) a 'frozen' or 'folded' state at low temperature and 

large stiffness OOjIIIl- We note that while the nexib le 
and semi-flexible F — T phase-diagrams are qualitatively 
similar [24| , the re-entrant behavior is suppressed because 
of stiffness and becomes less pronounced with increasing 
bending energy. In the CFE, the probability distribution 
of the end-to-end distance has "saw-tooth" like behav- 
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FIG. 5: Plot of (F) vs x for a semi-flexible chain with bending 
energy ef, = 0.3 at different T for N — 45 



ior corresponding to intermediate states during unfold- 
ing [24j ]. Therefore, it is important to study the effect of 
stiffness on force-extension curves in the CDE. The force 
extension curves shown in Fig. [5] have striking differences 
to the flexible ones. At low temperatures we see strong 
oscillations which vanish as the temperature is increased. 
Since the polymer chain has "frozen conformations" like 
(3— sheets (the zero-force limit of which describes zig-zag 
configurations inscribed in a square) [3Gj], it takes more 



force to unfold a layer. When about half a layer has been 
opened, the bending energy favors a complete stretching 
of the layers and hence the force decreases. This phe- 
nomenon allows us to probe a molecule like Titin which 
has similar /3-sheet structure [29j. 

To summarize, we have presented the exact solution 
of a model of long (finite) polymer chains, of direct rel- 
evance to recent experiments on the elastic properties of 
single biomolecules. The model takes into account several 
constraints imposed by the experimental setups: geomet- 
ric constraints, excluded volume effects, attraction be- 
tween chain segments, finite but large chain length (here 
up to N — 55). It permits one to choose the thermody- 
namic ensemble (CFE or CDE) dictated by the experi- 
mental protocol. The exact enumeration data permits us 
to access all parameter values, including biologically rele- 
vant low temperatures where previous studies have failed. 
Our results correctly reproduce several experimentally 
observed effects: multiple transitions during unfolding, 
"saw-tooth" like oscillations in the force-extension curve 
of semi-flexible chains and first-order transition into a 
"ball string" conformation. Finally, we have identified 
cross-over behavior that provides a unified treatment of 
both pulling and compressing forces in the CDE. 
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